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We propose to use the intrinsic two-level system (TLS) defect states found naturally in inte- 
grated optomechanical devices for exploring cavity QED-like phenomena with localized phonons. 
The Jaynes-Cummings-type interaction between TLS and mechanics can reach the strong coupling 
regime for existing nano-optomechanical systems, observable via clear signatures in the optomechan- 
ical output spectrum. These signatures persist even at finite temperature, and we derive an explicit 
expression for the temperature at which they vanish. Further, the ability to drive the defect with 
a microwave field allows for realization of phonon blockade, and the available controls are sufficient 
to deterministically prepare non-classical states of the mechanical resonator. 



Introduction. — Cavity optomechanics has en- 

abled the preparation of mechanical resonators in states 
of low phonon occupation via optomechanical (OM) side- 
band cooling [!-!)], and to observe their quantum co- 
herent interaction with light [ ]. OM systems have 
enabled displacement detection at or even below the 
standard quantum limit [10-12], thereby complementing 
other mechanics-based sensing applications [13, 14]. Fur- 
ther, they have been proposed for creating macroscopic 
quantum superpositions [15] as well as for applications 
in quantum information [16, 17]. In experiments carried 
out so far, the interaction between mechanical oscillator 
and cavity field is effectively linear, while one of the ma- 
jor challenges in the field is to reahze non-Hnearities at 
the single phonon level. The intrinsic OM radiation pres- 
sure non-linearity is predicted to enable the generation 
of non-classical states of light and mechanics [18, 19], 
provided that the single-photon coupling rate exceeds 
the mechanical frequency and the cavity decay rate. In 
multi-mode OM systems, the same non-linearity can be 
exploited more easily and it has been proposed to use 
it for enhanced readout [20] and quantum information 
processing [21]. 

Here we propose an alternative route to render the 
dynamics of the mechanical oscillator nonlinear at the 
single quantum level: using its natural coupling to in- 
trinsic structural two-level system (TLS) defects and 
thereby alleviating the need to functionalize the system 
[see Fig. 1(a)]. Ensembles of TLS defects were first stud- 
ied in the context of the anomalous and universal low 
temperature properties of glasses [22-20], where they 
arise from frustration. In experiments involving Joseph- 
son junctions, individual TLSs with transition energies 
distributed well into the GHz regime were observed and 
studied for their role in decoherence [27]. Nevertheless, 
their comparatively long coherence times, and their abil- 
ity to strongly couple to Josephson junctions via the elec- 
tric dipole moment have enabled a TLS quantum mem- 
ory [28]. In the same context, the influence of strain on 
TLSs has been probed recently [20]. However, in the 




FIG. 1. (a) Strain-coupling of a single two-level system (TLS) 
defect to an optomechanical system, (b) At low tempera- 
ture, the defect can be effectively described by two states in 
a double-well potential, where Ao is the tunnel splitting fre- 
quency and A the asymmetry frequency, (c) Schematic illus- 
tration of decay channels and couplings (see text). Resonator 
and TLS form a Jaynes-Cummings model (dashed box), ex- 
hibiting the characteristic anharmonic spectrum shown in (d) . 

OM setting, TLS ensembles have mainly been studied 
as a source of decoherence [30-32]. In this letter, we 
demonstrate theoretically that the coupling of an individ- 
ual TLS to a localized phonon mode of an OM system can 
be large enough to exceed the mechanical and TLS dissi- 
pation rates, and hence it paves the way for cavity QED- 
like experiments with single phonons. This interaction, 
as shown below to be described by a Jaynes-Cummings 
(JC) Hamiltonian [31, 33], induces single-phonon nonlin- 
earities that can be witnessed in the OM cavity output 
spectrum. Additionally driving the TLS with microwaves 
leads to phonon blockade, entailing a mechanical state 
with sub-Poissonian statistics; and more complex me- 
chanical states can be engineered using suitable proto- 
cols. Beyond this, OM experiments with single TLS may 
mature our yet incomplete understanding of such defects 
in glasses [■ U] . Our results apply equally to other classes 
of defect states - in particular, recent studies have con- 
sidered the strong coupling of GHz phonons to donor- 
acceptor impurity doped silicon [35, 3G]. 
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Strong coupling between resonator and TLS. — We con- 
sider optomechanical systems made of silica or silicon 
where a single optical cavity mode a couples to a co- 
localized deformational mode b of frequency ujrn- The 
Hamiltonian of the system including the intrinsic TLS is 
given by 



H — Horn + Hjc + -^TLS.p 



(1) 



Here, i/om = —h/S.^a'^a + hg{a + a^)(6 + W) describes 
the standard linearized OM coupling of rate g in a frame 
rotating at the frequency of the driving laser wz,, which 
is detuned from the bare cavity resonance by A/, [1, 2]. 
The remaining terms contain the interactions of the me- 
chanical mode with the TLS and the microwave drive of 
the TLS, as introduced below. 

Defects in low-temperature glasses are effectively de- 
scribed by TLSs with tunnel splitting /iAq and asymme- 
try energy fi.A [24], such that the eigenstates are split by 
At = v/A2 + A§ [see Fig. 1(b)]. Since A depends on the 
strain in the material, the TLS in the OM system cou- 
ples to the localized phonon mode b [24, 31]. The latter 
produces a zero-point strain fluctuation on the order of 

1/2 

5zpf = {fujJm/'^EVm) , where E is the Young's modulus 
of the material and the mechanical mode volume Vm is 
defined in analogy to optical cavity QED as [■ST] 



^ /T.,-(x)5,j(: 



ESij(xQ)Sij 



(2) 



Here, T^- and Sij are the stress and strain profiles, re- 
spectively, and repeated indices are summed over. Fur- 
ther, Xq denotes the point where Sij Sij becomes max- 
imal. Denoting by ct^ the Pauli matrices in the TLS 
eigenbasis, the interaction between resonator and TLS 
can be approximated by the JC Hamiltonian Hjc = 
(SAt/2)ct2 -I- fkOmb^b + hX{a+b + provided that 

A <C At « uJm [33]. Here, the TLS-phonon coupling A is 
given by [37] 



A 



h At 



huJm 

2EV„, 



1/2 



(3) 



where Dt is the deformation potential [38, 39] and we 
have neglected factors on the order of one due to the 
exact position and orientation of the TLS. 

In addition to strain, the TLS is also susceptible to 
classical electromagnetic fields [24, 4(J]. It responds to a 
coherent microwave drive of Rabi-frequency fJ^ accord- 
ing to i?TLS,M = ^f^/xe*'^"*'^- + fin*^e-"^>^*a+. On the 
other hand, a static electric field Eq causes a change in 
the asymmetry, SA = 2p ■ Eq/H, where p is the dipole 
moment of the TLS, and thereby changes the splitting by 
i5At = 2(A/At)p • Eo//i and the coupling by a smaller 
amount 5X = — 2A(A/Ay)p-Eo/fi. Also in the case of ar- 
tificial donor-acceptor based systems [36], similar tuning 
can be afforded by external electric and magnetic fields. 
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TABLE I. Finite-element simulations of mechanical mode pro- 
files for a radial breathing mode of a Si02 microsphere [11], 
a pinch mode of a Si spoke-anchored microdisc [42] , and a lo- 
calized breathing mode of a Si photonic crystal nanobeam 
[ ] (from top to bottom). Here Amax is calculated using 
Eq. (3) with Dt — 1.4eV [38], Ao = At, and well-known 
material properties [13]. Nt ~ O.SftAmaxV^m^ denotes the 
number of relevant TLSs per mode volume [37], i.e. those 
within a bandwidth A around ujm that have Aq/At > 0.7. 
P = 10*^ J"V"^ is the spectral density [:;i]. 



In Table I we display the results of full finite-element 
simulations of the mechanical modes of three different 
OM structures. In addition to cj„i, Vm and A, we also 
show the number Nt of TLSs that couple resonantly 
and appreciably to the mechanical mode (see caption for 
details). Clearly, a small mode volume entails a large 
coupling A oc l/VKn, but also reduces Nt oc \/Vm- 
While Nt ^ 1 is desirable in order not to couple to 
several TLS, a value of Nt <C 1 can be compensated 
by the above-mentioned tunability, which allows shifting 
an off-resonant TLS into resonance with the mechani- 
cal mode. In particular, already moderate electric fields 
[Eol - lO^V/m allow for shifts of 5At/2ti - 1 MHz, 
where we used |p| ~ 0.5 D [ ] and Aq/At ^ 0.9. The 
TLS-phonon couplings A estimated in Table I clearly ex- 
ceed the typical cryogenic mechanical linewidth 7m/27r 
10 kHz realized in recent experiments [ , ]. On the other 
hand, typical TLS relaxation rates 7T/27r have been mea- 
sured to be in the range 0.1 — 5 MHz for GHz frequen- 
cies [38, 44], and experiments suggest that they dominate 
over dephasing rates [44], which we thus ignore in this 
work. Using the bulk value 7t /27r 1 MHz at T ~ IK 
[ ], we obtain A/7T ^ 10 for the case of the nanobeam. 
The phonon density of states responsible for TLS relax- 
ation is strongly suppressed in the proposed structures as 
compared to the bulk, which leads us to consider these 
numbers a worst-case estimate. 

We conclude from the above discussion and the results 
in Table I that strong coupling to individual defects is 
feasible for suitably engineered OM systems. While the 
couplings in the microsphere and microdisc structures 
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are on the verge of the strong couphng regime, the most 
promising structure is the nanobeam. In the following, 
we derive the signatures of the TLS-resonator interac- 
tion in the OM output spectrum and show that the cou- 
pling enables quantum state control of the mechanical 
resonator in analogy to atomic cavity QED. 

Optomechanical output spectrum. — In the absence of 
the microwave drive (f2^ = 0), the system is only driven 
by thermal and vacuum fluctuations, as described by the 
master equation 

p = -l[H,p]+CcP + C,nP + CTP, (4) 
n 

with Lindblad terms CcP ~ K2?[a]p for the cavity, and 

CtP = 7t {nr + 1) 'D[a-]p + jTnT'D[cr+]p, (5) 

+ l)V[b]p + J,nnmV[b^]p, (6) 

for TLS and resonator, respectively [see also Fig. 1(c)]. 
Here, k is the cavity energy decay rate, nm,T are the Bose 
occupation numbers corresponding to an environment at 
temperature T and 'D[x\p = xpx^ — {x^xp — px^x) /2. We 
consider the case where the three systems are in reso- 
nance (— sa ~ At), and where the cavity adia- 
batically follows the resonator dynamics [n :s> g, 7T(n-T + 
l),7m("-m + !)]• In this regime, the OM coupling leads 
to cooling of the mechanical resonator and the TLS [.3.3] : 
For the case A ^ k of relevance in this work, the optically 
induced mechanical damping rate A^~'^ « / 1"^ can ex- 
ceed the corresponding heating rate v4(+) « g'^n/Aijj'^, as 
well as the rethermalization rate 'jrn'iT-m- In addition, the 
OM coupling transduces the resonator motion to the cav- 
ity output, and thereby enables the observation of the hy- 
bridized resonator- TLS subsystem by photodetection. In 
particular, we consider here the cavity output spectrum 
S{u}) = (K/27r)/^dr e-*("-"^)^(at(r)a(0))ss, where 
the angular brackets denote the average in the steady 
state of Eq. (4). 

In Figure 2 we display S{io) for frequencies around the 
blue sideband Wbiuc = i^l + as a function of tem- 
perature. The spectra show transitions between the lev- 
els of the JC model iJjc formed by resonator and TLS 
[dashed box in Fig. 1(c)]. In the case of exact resonance 
{ujm = At), the ground state energy of Hjc is zero and 
the excited states {ni) have energies ll)„± = nuJm ± ^\/n 
(n = 1,2,...), giving rise to the "JC ladder" shown in 
Fig. 1(d). For T — almost all population is in the low- 
est states of the ladder, such that only the transitions at 
uj — Wbiue i A between the first rung and the ground state 
are observed. As temperature increases, higher states get 
populated, so that transitions between higher rungs lo- 
cated at W = -f UJnaP, with iOnaP EE OJ^a - W(„_i)^ 

(a, /3 — ±), also become visible in the output spec- 
trum. Finally, above a certain cross-over temperature 
the spectrum consists of a single Lorentzian peak, since 
the upper and lower branches of the JC ladder resemble 
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FIG. 2. Fine-structure of the OM output spectrum S{ij) 
around the blue sideband at lumuc as a function of temperature 
and without a microwave drive field {Qfi = 0). The dashed 
line shows the cross-over temperature according to Eq. (8) 
and the lower panels show cross sections at the temperatures 
indicated by dotted lines in the top panels. Parameters are 

LOm. — 2-K-5 GHz, K = O.luJrn, A = 2 - 10~^CJm, 7m = 6- W~^LOm 

and (a) g = A, 7t = A/10, (b) g = A/10, 7t = A/30. 



two independent highly excited harmonic oscillators [4"i]. 
Note that below Tc, temperature actually helps to ob- 
serve transitions between higher rungs in the JC ladder. 

To estimate Tc, we adiabatically eliminate the cavity 
mode and calculate S{uj) using a secular approximation, 
which is valid as long as the individual spectral lines do 
not overlap [ ]. The resulting blue sideband of S{uj) can 
be approximated by a sum over Lorentzians li^{itJo,ja) — 



oo 



/ , / , ^na^ (w„a/3 , 7na/3 /2) : 



(7) 



where each term corresponds to a transition \na) — > 
\{n — l)/3) centered at uj = ujl + uJnaf)- Since A ^ k, 
the widths can be written as 7„a^ = 2(n — 1)7(^1 -I- 1) -I- 
2n7n-|-7T(2nm-|-l), for n > 2, and 7iq+ = 7(3n + l/2)-|- 
7T(2n™-|-l/2). Here, 7 = 7™-hA(")-A(+) is the efltective 
mechanical damping rate and n = {nmlm + A'^'^^)/^ the 
corresponding effective mean occupation, as known from 
standard OM cooling [4, "]. The weights in Eq. (7) are ex- 
pressed in general as = (2A(~)p„aS^^^)/(7r7„Q,;3), 
where Pna is the steady state population of the eigen- 
state \na). The JC matrix elements read B'^^p — 

(2n - 1 2a/3v^n(n - l))/4, for n > 2, and Bl^^ = 1/2. 

As temperature increases and population moves up the 
JC ladder, transitions between the zb-branches become ir- 
relevant, while the ones within each branch occur closer 
iouj — Wbiuc and thus contribute to the center of the spec- 
trum. We define the cross-over point as the temperature 
Tc where the separation of the two dominant Lorentzians 
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FIG. 3. Optically probed phonon blockade, (a) gi^\o) of 
the cavity as a function of drive frequency and temperature 
T for \Qfj,\ = O.IA. All other parameters as in Fig. 2(a). 
(b) ^c^' (0) of the cavity (solid) at temperatures indicated by 
dashed lines in panel (a), and corresponding (?m'(0) of the 
resonator (dashed). 



in Eq. (7) (at each side of uj 
i.e., where 2X{y/N - y/N - 



= '^biuc) equals their width, 
1) ~ with TV being 

the index n for which VF,^!|!:'^ is maximal. If we further 
restrict the parameters to a regime of experimental in- 
terest: < Jmn^ < A'--^ < 7t < A and > 1, 
where n'^ is the bath mean occupation at Tc, then the 
cross-over criterion yields [ ] : 



27Ty (1 + 3A(-)/7t)'^^' 



(8) 

Naturally, increases with increasing coupling A and 
decreases when increasing the TLS linewidth, but addi- 
tionally Eq. (8) shows that Tc can be increased by en- 
hancing the OM cooling rate A'-'K The above estimate 
for Tc agrees well with our numerical simulations, as can 
be seen from Fig. 2 (dashed lines) . For finite detuning 
\ojm — At I < A of TLS and mechanics the general pic- 
ture described above remains valid, although the spectra 
generally become asymmetric. 

Optically probed phonon blockade. — By driving the 
TLS with a weak coherent microwave field (|fi^| A) we 
can realize a phonon blockade in analogy to cavity QED 
[47] : if the JC system is cooled to its ground state and the 
drive is tuned to a transition |0) — > |1±), then the sub- 
sequent transition |1±) — >■ |2±) is suppressed provided 
that A ^ 7t, see Fig. 1(d). The resulting sub-Poissonian 
resonator statistics {gm{0) = {h^ b^ bb) / {h^ h)"^ < 1) per- 
sist at finite T if exp {~huj,n/kBT) ^ l^/^P/Tii+j which 
ensures that the thermal occupation is irrelevant com- 
pared to the one due to the drive. Since the cavity op- 
erator adiabatically follows the mechanics, i.e. a{t) « 
-2ig/K[b{t) - i{K/4uJmP{t)] 4- noise [-ii], we have that 



ffrn (0) ~ gc (0) = {a^ a'^ aa) / {a'^ a)^ to zeroth order in 
<C 1, such that the phonon blockade can be probed 

(2) 

optically. Figure 3(a) displays gc (0) of the cavity as a 
function of drive frequency and temperature. For low 
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FIG. 4. Preparation of non-classical resonator states based on 
Ref. [ ] . (a) Level scheme with available tools and (b) gen- 
eral sequence for preparation of an arbitrary state using OM 
cooling, strong TLS drive (|n^| ^ A), free evolution (~ A), 
strong TLS AC-Stark shift Sljt ~ - At| > A, and 

TLS-resonator detuning, (c) Fidelity = M\p{tf)\'i! m) as 
a function of temperature and for increasing M from top to 
bottom. Dashed curves: M = 1,2,3 and parameters as in 
Fig. 2(a). Solid curves: M = 1, 2, 3, 4, 7, 9 and parameters as 
before, except for 7t = A/30. The OM interaction is switched 
off after the initial cooling. 



T, one clearly observes the expected anti-bunching at 
uj^ = ujm ± A, while these features disappear for higher 
T due to thermal occupation of the JC ladder. As can 
be seen from the cuts presented in Fig. 3(b), the ^^^•'(O)- 
function of the cavity is an upper bound for the one of 
the resonator in the region of pronounced anti-bunching 
and therefore, an optical gc (0) < 1 indicates phonon 
blockade. 

Preparation of non- classical states. — The possibility 
to electrically tune and coherently drive the TLS, to- 
gether with the strong coherent coupling between TLS 
and resonator gives rise to the prospect of determinis- 
tically preparing quantum states by suitable protocols. 
As an example, we propose preparing the system close 
to its ground state by OM cooling and then generat- 
ing a non-classical resonator state by a scheme analo- 
gous to the one of Ref. [48]. The necessary sequence is 
illustrated in Figs. 4(a,b), and as an example we plot in 
Fig. 4(c) the resulting fidelities for preparing the state 
\^m) = (|0) + |M)) /\/2 {M = 1,2,...), in the presence 
of imperfections. Clearly, the fidelity decreases with in- 
creasing M and T. However, already for \/^t = 10 one 
can achieve J- > 0.95 for temperatures around lOOmK, 
which constitutes an exciting avenue for OM systems in 
existence today. 
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SI: COUPLING BETWEEN MECHANICAL 
MODE AND TWO-LEVEL SYSTEM 



For the modes of interest, we have free-surface boundary 
conditions, i.e., the stress must vanish on the surface dV, 



Starting from the classical description of a solid, we 
derive the Jaynes-Cummings type interaction of the lo- 
calized mechanical mode in an optomechanical system 
with an intrinsic two-level system (TLS) defect. 



Classical description of a linear and isotropic elastic 

solid 



The dynamical equation governing the motion of a lin- 
ear and isotropic elastic solid is given by [ ] 



(9) 



where the Einstein summation convetion is used and Ui — 
Mj(x, t) is the displacement field of the solid with respect 
to its equilibrium configuration. The quantities A and /i 
are the Lame constants, which completely characterize 
the elastic properties of the linear isotropic solid. They 
are related to the Young's modulus E and the Poisson's 
ratio v by 



A 



Ev 



E 



Equation (9) can be expressed more compactly as 



(11) 



where the elastic moduli tensor aijki satisfies the sym- 
metries a^jki = ajiki = aijik = akuj and is given by 



^ijki = ^SijSki + p{5ik5]i + SjkSii). 



(12) 



Also in terms of aijki, the generalized Hooke's law for a 
solid with linear response can be compactly written as 



Tij = OLijkiSkii 



(13) 



where Tij (x, t) is the applied stress field on the solid and 
Sij{yi,t) is the correponding strain field, which accounts 
for the degree of deformation of the solid, defined as 



Si 



- (d^uj + djUi) 



(14) 



Classical Lagrangian for a single mode 



The equation of motion (9) can be obtained from the 
Lagrangian 

L=\ j {pUiUi - aijkiSijSki) , (15) 

which, owing to the symmetries of aijki, can be simplified 
to 

L = ^ y d^x [puiu^ - a^jki{diUj){dkUi)\ . (16) 



(17) 



where rii is an outwards-pointing unit vector. 

Assuming the boundary condition (17), it is shown in 
Ref. [ ] that the differential operator p~^aijkidjdk on the 
right hand side of Eq. (11) is hermitian and therefore its 
eigenf unctions form a complete orthonormal set on the 
domain V . In our case, wc focus on a single, well-isolated 
and localized eigenmode of frequency aj,„, which can be 
described by the ansatz: 



Ui{x,t) = Ui(x)/3(i), 
such that Eqs. (9) and (16) reduce to 

= aijkidjdkUi + pool^Ui, 



L = 



1 



pu.iUiP'^ - aijki{diUj){dkUi)l3'^ 



(18) 

(19) 
(20) 

(21) 



Completing the total derivative in the second term of 
(10) Eq. (21), using Eq. (20), the symmetries of aijki and 
Hooke's law (13) for the spatial profile of the stress tensor 
Tij{x) = aijkidkUi{Ti), we are left with the Lagrangian, 



q2 



dPxUiUi 



2 .lav 



(22) 

This form makes explicit that for free-surface and/or 
fixed-surface boundary conditions, the system is simply 
described by an harmonic oscillator Lagrangian for /3: 



where the effective mass of the mode is defined as 



TOoff : 



(23) 



(24) 



Quantization of displacement and strain fields 

The quantization of the harmonic oscillator La- 
grangian (23) proceeds in the textbook manner and thus 
the Schrodinger picture operator /? can be expressed as 



2moffWm 



(25) 



where h is the destruction operator of the mechani- 
cal mode, satisfying [6,6^] ~ 1. The Hamiltonian op- 
erator describing the dynamics of b is simply Hm = 
hu^ {b^b +1/2). 
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Based on Eqs. (18) and (25), the quantized displace- 
ment field in the Heisenberg picture reads 



2TOoffW„ 



i,(x)(6e-*""-* + fo^e*""*), (26) 



such that the corresponding strain field obtained from 
Eq. (14) can be expressed as 



5,,(x,t) 



2moffa;,i 



5,j-(x) (6e-'"'"* + fe^e'""*) . (27) 



Here Sij{x) = {diUj + djUi)/2 is the strain profile of 
the mode of interest. We further introduce a normal- 
ized strain profile according to Sij(x) = Sij/Snorm, where 
Snorm = [^ij(p^o)Sij{xo)Y^^ , and xq is the point where 
SijSij becomes maximal. As a result, we have 

S,j{^,t) = 5,pfs„-(x) (6e-*""* + 6te*""*) , (28) 

where the zero-point strain fluctuation is given by 

/ hu!„ ^ ^^'^ 



S: 



zpf 



\2EV„, 



(29) 



In analogy to cavity QED [■>], it contains the mode 
volume Vm of the mechanical mode, which is defined 
as 1/™ = m^sLol/ESl,,^. Using Eqs. (13), (17), (20), 
(24) and the symmetries of aijki, we find the relation 
meffW^ = / (Px TijSki, such that the mode volume can 
be written as 



Vrr. 



/f„-(x)^,,(x) d^x 



(30) 



£;S'„(xo)5'jj(xo) 
This is the result displayed in Eq. (2) of the main text 



Hamiltonian for the interaction between defects and 
mechanical mode 

According to the theory of two-level system (TLS) de- 
fects in amorphous solids [ - ] , the Hamiltonian describ- 
ing a single isolated TLS is given by 



(31) 



where A and Aq arc the asymmetry frequency and the 
tunneling amplitude, respectively, and the ai are the 
Pauli matrices. The TLS couples to electric Eg and strain 
Sij fields via its asymmetry, i.e. A = A (Eg, Sij), and for 
weak fields we write 



we exploit the coupling to the electric field for driving the 
TLS with microwave radiation, but also to tune it by ap- 
plying a static field. 

The deformation potential term in Eq. (32) provides 
the mechanism for coupling the TLS to the localized me- 
chanical resonance of the optomechanical system intro- 
duced above. To derive the form of this interaction, we 
replace Eq. (32) with Eq = into Eq. (31) and use the 
quantum strain field (28) in the Schrodinger picture to 
obtain 

Ht ^^tiAa^-^hAoa^ + hX{b + b'^)a^. (33) 



Here the coupling A is given by 



1 



(34) 



where xy is the position of the TLS within the optome- 
chanical structure. We are interested in a situation where 
the eigenfrequency splitting Aj- = Aq + A^ of the un- 
coupled TLS is approximately equal to the resonator fre- 
quency in order to obtain a resonant coupling between 
the two. We further expect the coupling to be much 
smaller than the free frequencies, such that we concen- 
trate on the situation 



(35) 



It is thus convenient to change to the eigenbasis ("delo- 
calized picture" ) of the uncoupled TLS and to apply the 
rotating-wave approximation, yielding 



1 



(36) 



Here, the transformed operators read = (A/AT)cr^ — 
{Ao/AT)ax and = {Aq/At)^^ + (A/At)^^, and the 
general TLS-phonon coupling A is given by 



A = 



1 



Ao 
At 



/ huJri 



1/2 



AjSjj(xt). 



(37) 



Finally, neglecting factors of order one due to the exact 
position and orientation of the TLS, we can approximate 
Eq. (37) by 



A 



Dt 
IT 



Ao 
At 



2EV„, 



1/2 



(38) 



where Dt is the scalar deformation potential that has 
been measured in the literature (see, e.g., Ref. [9]). Equa- 
tion (38) is the expression for the TLS-phonon coupling 
displayed in Eq. (3) of the main text. 



A(Eo, S,,) « A + 2D,jS,, + 2 p • Eo, 



(32) 



where Dij = {l/2){dA/dSij) is the deformation poten- 
tial tensor and p = (l/2)(9A/9Eo) the instrinsic electic 
dipole moment of the TLS. As discussed in the main text. 



Estimating the number of relevant TLS 

The number of TLS with parameters around Aq and 
A is commonly modeled by the distribution cINt = 
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(fi.yP/Ao)dAo(iA, where V is the volume of interest and 
P the spectral density [G] . We are interested in the num- 
ber of TLS within a bandwidth A around the mechanical 
frequency ujm, since these can be considered resonant. 
Further, the TLS-phonon coupling in Eq. (38) should be 
appreciable, which means that we only count TLS with 
u = Aq/At above a certain lower cutoff uq. Setting 
the appropiate integration range and changing to more 
convenient variables, we thus evaluate 

Nt = VPh / du / dAT , (39) 

= VPhXn,,^ / du ^ , (40) 

where we used the fact that A depends on u according to 
A = Ainax", with Ainax = DrSzpi/h. Assuming uq = 0.7 
and identifying V with the mode volume Vm, the total 
number of relevant TLS is estimated to be 

Nt « /iA,„axKnP X 0.8. (41) 

This result is cited in the caption of Table I of the main 
text; it only moderately depends on the cutoff as long as 
uo > 0.5. 

SI: OPTOMECHANICAL OUTPUT SPECTRUM 

We provide details regarding the calculation of the op- 
tomechanical output spectrum in the absence of coherent 
microwave driving. In particular, we derive the expres- 
sion for the cross-over temperature Tc, at which the TLS 
signatures in the spectrum vanish. 

The Hamiltonian of the complete system including 
TLS, resonator, and cavity can be written as H — 
Horn + Hjc- Here, 

i/om = -hALa^a + hg{a + a^){h + 6^), (42) 

describes the laser-enhanced optomechanical coupling of 
rate g between the mechanical mode h and the cavity 
mode a, in a frame rotating at the frequency of the driv- 
ing laser [10, 11]. The detuning A^ of the laser from 
the bare cavity frequency lOc reads A^ = — lJc- In 
addition, the JC Hamiltonian 

Hjc = hurnb^b + ^HAtctz + h\{a+b + 6V_), (43) 

describes the interaction of the mechanical resonator with 
the TLS. Apart from coherent interactions we also ac- 
count for dissipative processes: The cavity is coupled to 
the electromagnetic vacuum giving rise to an energy de- 
cay rate k, whereas the mechanical oscillator and the 
TLS interact with indepedent heat baths at temperature 
T . The corresponding energy decay rates are denoted by 



7m and 7t, respectively. Again using A <C iOm ~ At, the 
dynamics of the total open system is described by the 
master equation 

p{t) = - '-[H,p\ + tzD[a]p 

~r '~im {rim + l)2?[6]/0 + 7„in„i2?[6^]p 

+ 7r(fiT + l)2'[a'-]p + 7tSt2'[o'+]/3, (44) 

where 'D[x\p = xpx^ — [x^xp — px^x)/2 denotes a Lind- 
bland term, and nm,T are the Bose occupation numbers 
of the baths, evaluated at the mechanical and TLS fre- 
quencies, respectively. 

Jaynes-Cummings eigenbasis 

The JC Hamiltonian iJjc can be exactly diagonal- 
ized [12]. Denoting by jn,!?) and \n,e) the states with 
n phonons and the TLS in its ground or excited state, 
respectively, the eigenstates \n±) of iJjc can be written 
as 

\na) ^Cna\n,g) + Sna\n-l,e), n > 1, a = ±, (45) 

and |0-|-) = jo, 5) for the ground state (|0-) = 0). The 
coefficients Cna and Sna are explicitly given by C„q — 

a [(ri„ -|- aSuj) /2U,n]^^'^ and Sna = [{^n — Oi5uj) /iVLn]^^"^ ^ 

where ri„ = (Jw^ + 4nA^)^/^ and the TLS-phonon de- 
tuning reads 5oj ~ ujm ~ A^. The JC Hamiltonian in its 
eigenbasis reads 

00 

Hjc = X! X! ^na\na){na\, (46) 

n—O Q— ± 

where the eigenfrequencies can be calculated to be 

w„a = rjw,„ -f -(ailn - (5w), n > 1, a = ±, (47) 

and uJo± = 0. Also in the JC eigenbasis, the operators b 
and cr_ can be decomposed as 

00 

b=Y.Yl BnapOnap, (48) 

n=l Q,/3=± 
00 

C^- = ^ ^ <7napOnal3, (49) 
n=l Q,;9=± 

where the transition operators between the eigenstates of 
adjacent rungs in the JC ladder are given by 

Onafi^\{.n^l)l3){na\. (50) 

The matrix elements CT„a/3 and Bnafi read 

Bnali ~ Ci^n-l)pCnaVri + '5'(ii-l)^S'„a\/'^ ~ I7 (51) 

(52) 

for n > 2. For n = 1, one has Bia+ = Ciq, CFia+ = Sia 
and the remaining coefficients vanish. 
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Effective master equation for mechanical mode and 
two-level system 

As discussed in the main text, we are interested in the 
situation where the optical cavity adiabatically follows 
the dynamics of resonator and TLS, giving rise to an 
optomechanical cooling effect. Based on the conditions 



(53) 



w„a^ being the difference between JC eigenfrequencies of 
adjacent rungs, uJnap = ^na — The dissipative 

dynamics of resonator and TLS is described by 

'^na/3n' a.' p' 

naj3 n'a.' ^' 

+ E E ^illn'o.'pMOiap,On'.'P')p, (56) 



nap n a , 



we thus derive an effective master equation for the re- 
duced density matrix ps of resonator and TLS by adia- 
batically eliminating the cavity mode. To this end, we 
transform the master equation (44) to an interaction pic- 
ture defined by the unitary U{t) = e-KH.,c-Ai,a''a)t 
express the result as 



(54) 



where p = pU and CcP = KT^lalp. Further, in the JC 
eigenbasis the optomechanical interaction reads 



Ant/5 = -igy ]Bnal3[ia + a^)iOnaP+Oia3),p], (55) 
na/3 

where the operators now carry a time-dependence accord- 
ing to a(t) = e'^-f-*a and Ona/^it) = e-'""'^^*Onap, with 



where 'D{A, B)p = ApB - {BAp - pBA)/2, and we have 
introduced the generalized transition rates 



+ 7t(^T + ^)(^nal3Crn'a'l3' , (57) 
^nal3n'a'l3' ~ ImnmBnafS Bn' a' 13' 

+ ^TnTCrnal3Crn'a'l3'- (58) 



Based on the separation of timescales in (53), we em- 
ploy standard projection operator techniques [l-i, 14] to 
derive an effective master equation for ps, valid up to 
order ~ Still within the interaction picture, we 

obtain 



Ps=CjcPs+g^^ E Brial3Bn'a'P''RG{R{uin'a'P')} (OnaPPsOl^'a'0' 



Oi^pdn'^'P'Ps+h.C. 



nct^ n'a' P' 



' 5^ E E BnapBn'a' l3''RG{R{-UJn'a' P')} (ol^a^Ps^"' a' f}' " OnapO^, ^, Ps + h.i 



nap n'a' p' 



E E BnapBn'a'P'lmlRiixJn'a'lS')} [On' a' P' PsO\afi + Ps^l^' a' ii'^naP ~ h.C. 



naf3 n' a.' P' 



nap n'a' , 



(59) 



where the response of the eliminated cavity enters via the To zeroth order in the small ratios Slu/k and X/k, we 



quantity 



R{oj) 



K/2+i{LU + Al) 

(/«/2)2 + (a; + Ai)2 



(60) 



can approximate R{±LUnap) ~ R{±u}m) in Eq. (59) and 
thus pull these factors out of the sums. By subsequently 
reverting to the operators b and cr_, we see that the ef- 
fective master equation assumes the same form as our 
original description of the dissipative JC model, i.e.. 



Effective Jaynes-Cummings description with 
modified parameters 



In addition to the conditions in (53), the setups of 
interest in this work can be assumed to satisfy 



<5a; < A < n. 



(61) 



Psit) = - -[Hjc,Ps]+lin+l)V[b]ps+jn:D[b^]ps 
n 

+ lT{nT + l)'D[a-]ps + jTnTT^[a+]ps. (62) 
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However, the mechanical parameters now incorporate the 
effect of the optical cavity: 



7 = 7™ + ^(-)-A(+), 



n 



7™ + A(-) - A(+) - 



(63) 
(64) 



where the optically induced cooling and heating rates are 
given as in Ref. [ ] by 



(65) 



Note also that we neglected a small mechanical frequency 
regularization of the order ^ juim <C in the JC 
Hamiltonian. 

The effective master equation (62) is a convenient 
starting point for the analytical and numerical analysis 
of the system, since it is simple in form in comparison 
to (59) and still describes all the physics we are inter- 
ested in. Furthermore, it remains valid when adding 
a weak coherent drive to the system, as described by 

ffTLs,M = nn^e'^^^'a^ + hnie-'^^'<7+, with « a. 

We therefore use it to obtain the results on the phonon 
blockade effect displayed in Fig. 3 of the main text. Here, 
a direct numerical treatment of the full problem is com- 
plicated by the fact that there is no interaction picture in 
which the full Liouvillian is time independent, such that 
(62) is extremely helpful. 

Calculation of optomechanical output spectrum 

The output spectrum of the optical cavity is defined 
as S{u) ^ (l/2n) C^dT (aLt(T)aout(0)).s [13], 

where aout(i) is the operator associated with the out- 
going field, and the angular brackets (...)ss denote the 
average in the steady state of the system. By using the 
input-output relation [13], the spectrum can be expressed 
in terms of the cavity mode operator: 



Re 



dre-'("-"^)"(a^(T)a(0))s 



(66) 



where we have omitted a contribution ~ 5{uj — ujl) due 
to the laser that drives the optomechanical system and 
causes the enhancement of the radiation pressure cou- 
pling. In the Quantum Langevin formalism, an adia- 
batic relation between the Heisenberg operators a{t) and 
Onai3 {t) can be derived by making use of the same sepa- 
ration of time scales already used for deriving the master 
equation (59). Omitting transients, we obtain 



a{t) 



igBnapOnapit) 

^ k/2 - i(a;„Q^ + 

nap ^ 
t 



^ k/2 -f i(a;„Q/3 - Al) 

neap ' ^ ^ 



(67) 



Replacing this expression in Eq. (66), we see that we need 
two-time correlations functions between different transi- 
tions operators, as for example {0\^p{T)On,'a' pi{^))ss- 

For strong coupling and low enough temperatures 
the spectrum consists of non-overlapping peaks, each 
of which is associated with a transition between JC 
eigenstates. This situation is well described by the sec- 
ular approximation [lo, 1(>], in which cross-terms be- 
tween different transition operators are neglected based 

on the assumption \LOnaii - ^n'a'P'\ > W^^iin' a' Si-V for 
{na(3) {n'a' (5'). Applying this approximation to (59) 
yields a result similar to the one obtained in Ref. [IG]: 

Ps = -i[Hjc,Ps\ 

+ E t'^-^] + ^n>PL(3]Ps} , (68) 

nQ/3 

where the transition rates read 

= l{n + 1)BL/5 + lT{nT + l)aL^, (69) 



^(+) 

" na/3 



(70) 



Note that Eq. (68) can also be obtained by applying the 
secular approximation directly to Eq. (62). 

The simplicity of (68) carries over to the equations of 
motion for the expectation values: 

|(OnoMi)> = - {J-^+^^n»p) {OnMt)), (71) 

where the effective decay rate for the transition with fre- 
quency ujnais reads 



7nQ/3 = 



(+) 



(72) 



By means of the quantum regression theorem [13, 14], 
this result enables us to calculate the two-time corre- 
lation functions needed to evaluate (66). As a result, 
the spectrum can be written as a sum over Lorentzians 
L(wo,7o) = 7o/(7o + (t^ - - t^o)^): 

S{uj) = 5i.od(w) + 5biuc(w) (73) 

+ E ^-1 ('^-^' ' (74) 

where we have separated the contributions of the blue 
and red mechanical sidebands of the spectrum, introduc- 
ing the weights: 



red 



w 



nap 



(75) 
(76) 
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In Eqs. (75)-(76), p„a = {\na){na\)ss denotes the steady 
state mean occupation of the JC eigenstate \na). The 
blue sideband term of Eq. (74) is the resuU quoted in 
Eq. (8) of the main text. 

Using Eq. (68) we find that the occupations Pna satisfy 
a classical master equation, whose steady-state is deter- 
mined by the recursion: 



+ Erit>(«-i)M = 0, n>0, (77) 



A" 



and where for compactness we defined F^^^ =0, Fq^^ = 

0, po- = and P(-i)a — 0. Since F^^^^ — F^^'^, it is easy 
to see that detailed balance solves Eq. (77) for n = 0: 



Pi+=Pi_ = (f1+VM+^+,Po+. 



(78) 



Also noticing that the rates in Eqs. (69)-(70) satisfy 

ritja = ^^nfia *^nd that the sums Y,,^ rla^^ are indepedent 
of a, it is straightforward to show that Eq. (77) implies 

Vn > 1. (79) 



Pn- 



Pn+, 



As a result, Eq. (77) simplifies and can be readily solved 
by detailed balance, yielding [ I ''] 



P{n+1) + 
Pn+ 



"/n{2n -|- 1) -|- jTi^T 
7(n + l){2n + 1) + jrinr + 1) 



(80) 



for n > 0. In order to fix the ground state occupation 
Po+ one must impose the normalization condition po+ + 
2 Pn+ = 1 and by this the complete steady state is 

obtained. 



Estimation of cross-over temperature 

Our final goal is to estimate the cross-over tempera- 
ture Tc, at which the individual peaks in the optome- 
chanical output spectrum S'biuc 

(w) visible for T < 
merge into a single resonance at w = t^biuo ^ + ^m- 
Within the secular approximation and A ^ k, as in the 
previous subsection, the spectrum is always symmetric 
around Wbiue and consists of many overlapping peaks as 
Tc is aproached. In particular, only transitions within 
the upper and lower branches of the JC ladder con- 
tribute, and we focus on the two most dominant peaks 
on either side of Wbiuc, corresponding to the transitions 
|n— ) — > |(rt— 1)— ) and \n+) — )• |(n— 1)+), with maximum 
weights iy^++ = W^^^il. Here N denotes the index n 
for which the weights are maximal at given T. Then, the 
cross-over approximately occurs when these two highest 
Lorentzian peaks appreciably overlap, i.e. when the dis- 
tance between their centers is approximately the same as 



the sum of their half widths: 

^ (7W++ 



■ 7JV- 



/2. 



(81) 



In the case of exact TLS-phonon resonance Suj — and 

since a;jv++ — ujn = 2a;jv++ and — Jn++, the 

cross-over condition (81) can be simplified to 



2X{VN - VN-l) « 7Ar4 



(82) 



In order to find the index N = N{T) as a function of 
temperature, one needs to maximize W^,^^+ in Eq. (76) 
as a function of n for given T. Since a simple analytical 
expression for the ratio P{n+i)+/Pn+ is given in Eq.(80), 
it is convenient to define /i„ = such 
that N can be obtained from the condition /ijv ~ 1, i-e.. 



P(N+1) 



R2 



1N++ 



Pn^ 



^N++ 7(Af+l)++ 



1. 



(83) 



Therefore, by using Eqs. (51), (72) and (80), one can in 
principle solve the above equation for N. However, before 
doing so, we can make further simplifications by restrict- 
ing the parameters to the regime of experimental inter- 
est. Denoting by the mechanical Bose occupation of 
the bath at the cross-over temperature we write the 
conditions: 



« A(- 



< 



7t, 



> 1. 



(84) 
(85) 
(86) 



Here, (84) are the conditions for optomechanical ground 
state cooling, while (85) states that the TLS decay rate 
is the largest decoherence rate. In addition, we ex- 
pect that for strong coupling A 3> 7t, the cross-over 
happens at temperatures corresponding to large mean 
occupations of the bath, which is expressed by (86). 
The above conditions imply that 7 ~ A^"' < and 
n w (7mn^)/^(^' <C 1 ^ As a consequence, we are 
able to perform a Taylor expansion of Eq. (83) in powers 
of {NA'^^'>)/{jTnm) ^ 1, which to zeroth order gives a 
quadratic equation for TV: 



c{c-lfN^- 



O 



fNA^- 



= 0. 



(87) 



Here c is a dimensionless constant of the order one, given 
by 



c= 1 



1 



2v4(-) 

"m7T 



O 







2 


1 ■ 




{ 


.7Tri^. 









With the help of (88) we can easily solve (87) for N, 
yielding 



N = 



1 



1 + 2A(-) /7 
l + 2A(-)/75 



1 



> 1. 



(89) 
(90) 
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Finally, we can plug our estimate (90) in the cross-over 
condition (82) and solve for n^. By applying the approx- 
imations (84)-(86) to 7JV+-I- given in (72) and assuming 
from Eq. (90) that iV 1, we are left with the equation, 



A(-)/7j 



1 + 2A(-)/7T 



A 



2VN 



If we now square Eq. (91) and use expression (90), we 
obtain the following analytical estimate for nj^: 



_X_Y' (1 + 2A(-V7t) 
27Ty (i + 3A(-)/7t)'/'' 



(92) 



This result consistenly satisfies our initial assumption, 
^ 1, provided the strong coupling condition, A 7t, 
is fulfilled. By approximating the Bose distribution for 
large occupations nj^ ^ 1, the corresponding cross-over 
temperature reads 



nt,-, 



(1 + 2A(-)/7t) 

27T, 



(1 + 3A(-)/7t)'^'' 
(93) 

which is the result quoted in Eq. (9) of the main text. 
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